Bosonization in arbitrary dimensions 
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Using methods of functional integration, and performing simple Gaussian integrals, I show that an 
interacting system of electrons can be bosonized in arbitrary dimensions, in terms of the electrostatic 
potential which mediates the interaction between them. Working with the bosonic field, the sytem 
is shown to exhibit localized structures reminiscent of striping in the cuprates. 
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Attempts at bosonizing a gas of interacting fermions 
date back to Tompnogald, in the context of superconduc- 
tivity. Luttingera followed with an exact solution of a 
model Hamiltonian in one dimension. Licb and Mattis 
perfected Luttinger's solution shortly thereafter. Luther 
and HaldaneQ published seminal papers which extended 
the older one-dimensional methods to higher dimen- 
sions. Their ideas iave been developed in recent years 
by Houahton ct alcl, Feldman et alQ, Castro-Neto and 
FradkinH, and Frohlich and GotschmannB. Other theo- 
retical methods involving techniques to bosonize fermions 
in one ot_twp dimensions have also appeared in the 
literatuxeEnia. These ideas almost exclusively entail the 
use of fermionic creation and annihilation operators, and 
as such are difficult to apply in greater than one spatial 
dimension. 

One of the chief motivations for bosonizing a system 
of interacting electrons is to investigate the possibility of 
these electrons forming bound states which may not be 
accessible via the usual perturbative methods. Recently, 
a Luttinger liquid has been invoked by AndersonE£l in 
an effort to understand high temperature superconduc- 
tivity (HTSC). We will provide in this paper a simpler, 
exact alternative to the Luther-Haldane formalism. The 
main focus of this paper will be to establish a straightfor- 
ward way to replace a system of interacting electrons with 
an equivalent system of interacting bosons using func- 
tional integral methods. In this paper it will be shown 
that bosonization of a certain kind can always be applied 
exactly to interacting fermionic systems in arbitrary di- 
mensions. The ideas established in this paper may lead 
to insight into the mechanism of high temperature su- 
perconductivity, especially as the detailed physics of the 
HTSC compounds is included into the problem. This 
paper should also be applicable to strongly correlated 
electron systems in general. It should also be relevant 
to the study of striping in the cuprates. In fact we will 
provide later in the paper a qualitative model for this 
phenomenon of striping. 

Bosonization has attracted attention in work on rela- 
tivistic field theories as well. In this paper I shall carry 
out in the non-relativistic regime appropriate to con- 
densed matter physics, the program of bosonization de- 
veloped recently by Bannerjee in relativistic Quantum 



Field Theory (QFT)00. Bannerjee used the method of 
functional integration to eliminate matter fields in the 
standard expression for the action, leaving just bosonic 
fields for consideration. The focus in this paper will be 
on phenomena which emerge when these methods are ap- 
plied to issues in condensed matter physics for d > 1. 

Our technique is related closely to .-the work of 
Fogedbyll3, Lee and CherO and KopietzEl The focus 
in these papers has been the computation of the elec- 
tronic Green's function, in an effort to understand the 
origin of non- Fermi- liquid behavior. I shall instead con- 
centrate on uncovering physics which lies in the effective 
action for the bosonic fields that replace the fermionic 
fields. As indicated above, I obtain a tentative model for 
striping in the cuprates. This paper shows how to elim- 
inate the approximations invoked by Kopietz to carry 
out successfully the program of bosonization in d > 1. 
Last but not least, the techniques utlizied here will be 
seen to be more economical than in previous papers, as a 
result of which the final, novel results are more transpar- 
ent. Rather than diminishing previous results, it is hoped 
that the perspectives of this paper can be combined with 
existing techniques to shed more light on HTSC. 

Without further ado, I note that the usual partition 
function for a gas of interacting electrons (in atomic 
units; h = c = e = m = 1) is given by 

T>ip*T>il> exp iS 

dt J d d xvb*(x, t) (id t - H ) t) 
- / dt ( d d xip*(x,t) [ d d x'\ip{£ ,t)\ 2 \x-x'\- 1 i;(x,t) 
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where ip and are two- vectors describing fermion 
fields having a spin 1/2. Hq can be modified with an ex- 
tra term describing the interaction of the electrons with 
an underlying lattice. This partition function can be de- 
rived by introducing the scalar electrostatic potential: 

Q ~ / V4)*V4>V4> exp iS' 
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S' = dt d d 



i>* (id t -4>)ip - ip*H tp 



8tt 



(2) 



where an irrelevant normalization constant has been 
ignored. 

Integrating out the electrostatic potential using Gaus- 
sian integration in || leads directly to |l| On the other 
hand, if the matter fields ip and tp* are integrated out 
in ||, the partition function involves just the electrostatic 
field <fi. In general, the electromagnetic vector potential 
A should also be included. Traditionally this is not done 
in condensed matter physics on grounds that it is the 
Coulomb interaction which dominates the physics. Elim- 
ination of the matter fields then leads to: 



Q ~ / Vcfrexpi dt / d d x { - 
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(det [_i(i$-0 + M-H)]) _ 



(3) 



where fi signifies a constraint used to guarantee num- 
ber conservation, since the integration over the matter 
fields includes cases involving broken symmetry. Note 
that this constraint is formally equivalent to an inter- 
action of the electrons with a uniform background field. 
The way /x has been defined, [i > implies an interaction 
with a uniform positive background. We see here an anal- 
ogy with the Wigner-Seitz model. An increasing number 
density of electrons is signalled by an increasing /i. The 
appearance of the scalar field cj> in ^| implies that we are 
dealing with a spin-zero object, but it should be remem- 
bered that <j) is an approximation to the photonic field, 
which carries a unit spin. Given that the gauge degree of 
freedom represented by <fi can be absorbed into the elec- 
tronic wavefunction, we see that there is a certain sim- 
ilarity between our approach and that of Haldaneu who 
performed a decomposition of the electronic wavefunction 
in terms of its phase. I will now evaluate the determinant 
in |[ which contains non-trivial physics. This evaluation 
is more convenient in Euclidean space, which entails set- 
ting t — > —it. I then need the determinant of the opera- 
tor (d t + <j) — n + Ho). Upon considering an associated 
heat diffusion equation viz., —d T g(x — x',t — t',T) = 
(dt + 4> — ft + Hq) g (ip,?, t — t 1 , t) the determinant is 
obtained as exp — C'(0)El, where 
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d 3 xg(6,0,r) (4) 



Here attention has been restricted to three dimensions 
for ease of presentation. Our method of evaluating de- 
terminants of operators is fairly well-establishectj. The 
basic idea is to perform the evaluation formally with a 
<j) which is constant, and then to use the form obtained 
thusly to write down the complete expression for the ef- 
fective potential. The manipulations of this equation are 
fairly straightforward, and the only care one needs to ex- 
ercise is in taking Laplace transforms, considering only 



those frequency domains which ensure convergence of the 
overall expression. The effective potential which arises 
from the determinant is given by: 



K//(^-m) = (I0-mI) 5/2 x 

[*(-5/2)-r(-5/2)ln(|0 



Ml)] (5) 



where in g is a digamma function and should not 
be confused the fermionic field used earlier. Given that 
we have here a nonrelativistic formulation of electron dy- 
namics, it is seen that temporal derivatives do not ap- 
pear in the bosonized version of the partition function. 
(j) refers to the (time-independent) electrostatic poten- 
tial. As such the partition function involving the field (j) 
can be redefined with a delta-function, in much the same 
way as gauge constraints are imposed on path-integral 
formulations. On classical grounds only > needs to 
be considered, since we have a system of interacting like 
charges. More generally, if we are considering quantum 
effects, negative values of <f> could occur. It can be shown 
graphically that for sufficiently small values of /x, the ef- 
fective potential displays a single minimum for > 0. 
There is of course a second minimum for negative val- 
ues of the field. I will revisit this double- well form of 
V e ff shortly. The compact form for V e ff was obtained 
relatively economically. Besides the ease of derivation, 
the simple form for V e f / displays a double- well structure 
which will shortly yield some novel insights into the in- 
teracting electron system we started out with. 

If the potential is expanded around the positive mini- 
mum which occurs for sufficiently low values of the chem- 
ical potential, then the Euler-Lagrange equation associ- 
ated with the effective Lagrangian can be written as: 



V 2 - A- 
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This is the analog of Laplace's equation for a system 
interacting via a screened Coulomb interaction. The right 
hand side of |^ signifies that the field configuration is such 
that there are no free charges in the system. I interpret 
this as an insulating phase! The reason for restricting 
attention to positive values of the field is that the Euler- 
Lagrange equation represents the classical limit. The ba- 
sic assumption behind this argument is that the system 
settles into a state desribed by </> m j„ and harmonic per- 
turbations around it. 

More complicated situations than the one described by 
can also occur. The previous discussion regarding the 
insulating regime can be expanded by considering the fol- 
lowing double- well representation for the effective poten- 
tiaLviz., V e ff = y o -(a/2)(0-/i) 2 + (/3/4)(0-ix) 4 . The 
two minima for this function occur at (f> m i n = /x± \J (xj ' (3. 
From this it follows that the Euler equation is: 



I 







(7) 



2 



where <j> = <f> — fj,, and £ = \J Ana x. Note that for /x 
sufficiently small, there is only one mimimum for > 0. 
In one dimension, as x — > 0, w (l/2)a; 2 (— (/3/a)/x 3 ), 
while for a; — * ±oo, <f> ~ 4> m in- This allows us to con- 
nect neighboring insulating regions via a domain wall. 
Going to a cylindrical geometry, it is possible to show 
approximately that as the cylindrical co-ordinate p — ► 0, 
cf> — fx ~ p s , with s=0,l,2... While for p — > oo, — (J, ~ 
0mm [l — •s 2 /(2/0 2 )] . By the arguments given above, the 
domain wall between the defect and the surrounding in- 
sulating phase has free charges flowing through it. An ar- 
ray of such localized defects are reminiscent of the strip- 
ing which occurs in the cuprates. The model we just 
described is-, complementary to the approach of Castro- 
Neto et allij. In spherical geometry, as the radial co- 
ordinate r — > 0, 4> — fi ~ r l , with 1=0,1,2,... For dis- 
tances far from the center of such spherulitic defects, 
(j) — [i ~ (f> m in [l — 1(1 + l)/(2r 2 )] . This approximate 
treatment in non-planar geometry can be trivially gener- 
alized using the unabridged form of V e f / . 

As \i increases, both minima in V e ff occur for pos- 
itive values of the field tj>. As indicated below |^, the 
appearance of fi is formally equivalent to introducing a 
Wigner-Seitz positive uniform background model. We 
would therefore like to argue that for sufficiently large 
values of (i, the system behaves as a metal. We cur- 
rently cannot describe this transition from an insulat- 
ing to a metallic phase. The only remark we can make 
is the following: Given that classically the system can 
settle into one of two energy minima, it must be pos- 
sible to have spatial domains in which the system has 
4> = ±(f> m i n , with domain walls permitting a transition 
from one state to the other. For example, if we assume 
that the three-dimensional system exhibits variation in 
only one direction, then the domain wall can be described 
by <p — H + 4>min tanh[x\f2~Tva\. This domain wall has free 
charges flowing through it, i.e., it satisfies Q, but with a 
source term on the right hand side. But this is similar to 
the defects discussed above, in the insulating phase. It 
is only if there are a large number of these domain walls, 
and they form a connected network, will the system con- 
duct appreciably. 

We are currently in the process of computing the elec- 



tronic Green's function, given that the system of inter- 
acting electrons can be bosonized. 

I would like to acknowledge electronic correspondence 
with P.W. Anderson and P. Kopietz rectifying gaps in 
my knowledge of previous literature. 
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